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Abstract 



We present the optimal relay-subset selection and transmission-time for a decode-and-forward, half-duplex 
cooperative network of arbitrary size. The resource allocation is obtained by maximizing over the rates obtained 
for each possible subset of active relays, and the unique time allocation for each set can be obtained by solving a 
linear system of equations. We also present a simple recursive algorithm for the optimization problem which reduces 
' the computational load of finding the required matrix inverses, and reduces the number of required iterations. Our 

O , results, in terms of outage rate, confirm the benefit of adding potential relays to a small network and the diminishing 

marginal returns for a larger network. We also show that optimizing over the channel resources ensures that more 
relays are active over a larger SNR range, and that linear network constellations significantly outperform grid 
^ ' constellations. Through simulations, the optimization is shown to be robust to node numbering. 



I. Introduction 

Cooperation has become a popular technique to implement diversity in the absence of multiple antennas 
at receiving and transmitting nodes [l]-[3]. In this context, resource allocation in cooperative networks 



^ has recently become an active research area, and has been investigated under many scenarios and metrics. 
In this paper, we address the problem of resource allocation, in terms of channel resources (time or 
bandwidth), in multi -relay networks with arbitrary connections. We describe the contributions of the 
paper in detail after a brief review of the pertinent literature. 

For the single-relay case, several works have dealt with various aspects of resource allocation, in terms of 
power and/or bandwidth and time. Yao et al. determine the optimal power and time allocation for relayed 
transmissions specifically in the low-power regime [4]. Larsson and Cao present various strategies for 
allocating power and channel resources under energy constraints [5]. For the channel resource allocation 
problem, however, the authors consider selection combining only and do not address the scenario of joint 
decoding of the source and relay signals. The works in [6]-[8] address the problem of power and channel 
resource allocation under sum average power constraints. Optimal time and bandwidth allocation using 
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instantaneous and average channel conditions is obtained using power control in [9]. Channel resource 
allocation under fixed power is developed in [10] 

In networks with multiple relays, the available literature can be classified into two groups: networks 
where relays do not communicate with one another (parallel-relay networks), and networks without 
restrictions on relay communication (arbitrarily-connected networks). Resource allocation for the for- 
mer has been addressed in [11]-[14]. Ibrahimi and Liang develop the optimal power allocation for 
a multi-relay cooperative OFDM A amplify-and-forward (AF) system [11]. By maximizing the channel 
mutual information, Anghel et al. find the optimal power allocation for multiple parallel relays in AF 
networks [12], [13]. A more general solution is given in [14] where the authors give the optimal power 
and channel resource allocation for a parallel-relay network with individual power constraints on the 
nodes. 

To the best of our knowledge, the problem of channel resource allocation for arbitrarily-connected 
networks and dedicated multiple access has not been addressed in the literature. In general, works in 
the area of multi-relay systems with arbitrary links generally neglect the bandwidth penalty arising from 
multiple hops by assuming either full-duplex nodes, a bandwidth-unconstrained system, or the availability 
of channel phase information at the transmitter [15]-[27]. 

These assumptions, however, are not realistic for practical wireless networks, where nodes are likely 
to be half-duplex, phase information is very difficult to obtain at the transmitter, and bandwidth is a 
scarce resource. To fill this void, in this paper we investigate the problem of resource allocation in a 
bandwidth-constrained, cooperative, decode-and-forward (DF), wireless network, and consider the most 
general setting where multiple relays can transmit can cooperate with each other in transmitting information 
between source and destination. In this setting, we address the joint problem of optimal selection of a 
relaying subset and allocation of time resources to the selected relays. The resource allocation is framed 
in the context of mesh networks of relatively simple and inexpensive nodes. We concentrate on resource 
allocation in terms of transmission time only, removing power allocation from the optimization; we further 
simply the problem by considering orthogonal transmissions. This is motivated by the need to reduce 
complexity, allowing for nodes which can implement the resource allocation simply by switching on 
and off. For a system without power allocation, a solution to this problem provides an upper bound on 
cooperative performance in multi-relay network where dedicated channels are assigned for each source 
transmission. 
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To the best of our knowledge, no other work provides a solution to time-allocation for an arbitrarily 
connected cooperative network. The solution can be interpreted as a generalization of the opportunistic 
protocol presented by Gunduz and Erkip, where the relay is active only when it increases the outage 
rate [6]. In terms of the resource allocation solution, it is also a generalization of the solution in [10], where 
channel resource allocation is determined under fixed power for a three-node DF network. The solution 
can also be interpreted as a generalization of node selection [3], [28]-[30] under relaxed transmission 
constraints, where transmission can occur on multiple time-slots and relays can communicate with one 
another. 

This paper is structured as follows. Section |ll] describes the system model. In Section Un] and |IV] 
we develop the proposed resource allocation scheme and present a significantly simplified recursive 
implementation. Simulation results are presented in Section |V] and concluding remarks are presented 
in Section |VIl 

II. System Model 

We consider a mesh network of static nodes comprising a source and destination node and potential 
relays. The inter-node channel powers are denoted as where i and j represent the source node s, 

relay nodes r^, k = 1 . . . N, or the destination node d. They are assumed independent of each other and 
are modelled as flat, slowly-fading and exponential with parameter A. A is inversely proportional to the 
average channel power and is a function of inter-node distance, dij, through the path loss exponent pa, e.g., 
^/Kd = i^/d^sd)^ ^I^Tkd = (l/'^rftd)- ^'^^ includc shadowing into the fading model, although 

this can easily be incorporated on an instantaneous basis. Because the nodes are static, the channels are 
assumed to change very slowly with time; we thus assume knowledge of all channel gains (although not 
channel phases) at a centralizing unit. This knowledge is essential to our resource allocation scheme. 

With the aim of designing simple and cheap nodes, we assume half-duplex channels and orthogonal 
transmissions, which greatly simplifies receiver design. The relays are assumed to be numbered in some 
convenient order such that relay Vj transmits after if j > i. For example, the relays may be in a 
linear constellation as shown in Figure \\\ We also assume the DF cooperation strategy with independent 
codebooks, which allow for the optimization of system resources (see [31] for an overview on current 
coding methods for nodes using DF). Note that repetition coding does not allow for this resource allocation. 

With these assumptions, the cooperation framework for the A^-relay fully-connected network is as 
follows. The half-duplex constraint precludes the relays from transmitting and receiving simultaneously 
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on the same channel, and the unavailabiUty of forward-channel phase information at the nodes precludes 
the nodes from simultaneous transmissions. The transmission between the source and destination is thus 
divided into + 1 time-slots, of duration to, h, . . ., In, with to + + • • ■ + ^at = 1. In the first time-slot, 
the source transmits its information to all the nodes. The first relay, ri, decodes this information and the 
remaining N relays and the destination store the information for future processing. In the second slot, 
of duration ti, the first relay re-transmits the information using an independent codebook, the second 
relay decodes the information from the first relay and the source, and the remaining N — 1 relays and 
the destination store the information for further processing. In general, each relay decodes information 
from the source and from the previous relays ri . . . i up to and including time-slot tk-i- This process 
continues until all relays have transmitted and the destination attempts to decode the information. 

Assuming that each node uses power P and W Hz per transmission (noting that although each node 
transmits for a different amount of time, the symbol durations and thus the corresponding bandwidth used 
by each node is the same), the signal to noise ratio (SNR) at node j resulting from transmission from 
node i can be written as SNRjj = jv^l^y P' where A^o is the noise spectral density. In the rest of the 
paper, we use the short-hand notation Lij to denote log2(l -I- SNR^j), the capacity of the corresponding 
channel. 

III. Optimal Resource Allocation and Relay Selection 

In this section, we solve the joint problem of resource allocation and relay selection for the network 
discussed above. Essentially, we give the optimum values of U, i — . . . N, such that the achievable rate 
between source and destination is maximized. We begin here with a fully-connected network, where each 
node is linked to all other nodes through a non-zero channel. 

A. Fully Connected Network 

Consider a source-destination pair communicating with the help of N relays. Assuming that each relay 
is active, the mutual information at each relay and destination can be written as 

-^i(^o) = toLsn, (1) 

Ik{to,ti,t2, ■ ■ ■ ,tk-l) — to-^srfc + + ■ ■ ■ + ^fe-l-^r-fc_ir-fc, (2) 

Id {to, ti, t2, ■ ■ ■ ,tk, ■ ■ ■ , tisf-l, tN) — toLsd + tiLr^d + . . . + tk-lLrf._-i^d + • • • + tiSfLrNd, 



(3) 
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where 7^ and Id denote the mutual information at relay and the destination, respectively. 

With all N relays cooperating, the maximum achievable rate under orthogonal transmissions is the 
minimum of the mutual information obtained at each individual relay node: 

Rn = max min {/i(to), ^2(^0, ^i), • • • , -^ik-i(^o, ■ ■ ■ , -^A;(^o, ■ ■ ■ , ^fe-i), -^A;+i(^o, • • • , ^it), 

to,.--,tjv 

. . . , lN{tQ, ■ ■ ■ ,tk, ■ ■ ■ , tN-l), loito, . . . ,tk, . . . , ^iV-l, t^)}, (4) 

such that ti> 0, Wi, 

to + ti + ..An <1. 

The above expression is a straightforward generalization of the cut-set bound for the single-relay network. 
This generalization maintains orthogonal transmissions for each relay, a model which represents practical 
networks with simple nodes that cannot implement complex interference cancelation. We use this model as 
the basis of the optimization in the rest of this paper. We note, however, that because each relay transmits 
using an orthogonal channel, Rn is clearly not the channel capacity. For literature on the channel capacity 
of arbitrarily-connected networks, we direct the reader to [32]-[37] for full-duplex relays, and [38] for 
half-duplex relays. 

For reasons that will soon become clear, let us consider the case with relay removed from the 
network. The maximum achievable rate Rn-i becomes 

Rn-1^ max min {/i(to) • • • , 4-i(to, • • • , ^-2), 4+i(to, • • • , ^fc-i), (5) 

to,---,tk-l,tk+l,---tN 

. . . , loito, ■ ■ ■ , tk-1, tk+1, . . . , ijv)}, 
such that U >0, Vi, 

to + . . . tk-i + ife+i + . . .In 

Removing relay is thus equivalent to removing tk and 1^ from the optimization. [We use the subscript in 
R%-i to denote the maximum number of potentially active relays, and the superscript to denote the relay 
removed]. The maximum rate at which the source can transmit to the destination can thus be written as 
the maximum of the rate obtained by using all relays, and the rate obtained by successively removing 
each relay: 

Rt — max{i?jv, -Riv-i, Rn-i^ ■ ■ ■ ■> -^jv-i}- (6) 
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If Rt = Rn^i^ the maximum rate can be obtained by iterating through (H)) and ([5]), successively removing 
a relay each step. Note that obtaining R%^i includes the cases where two or more relays are removed. 
In theory, therefore, all 2^ possible cases must be checked. 

Let (to)^i) • • • ^'t*N) denote the resource allocation that solves the optimization problem. We begin an 
outline of the solution to the optimization problem in dH), (|5]) and ^ with the following proposition. 

Proposition 1: With a maximum number of potential relays A^, the maximum achievable rate Rt = Rn 
only if tl ^ 0, VA;. Otherwise, if tl = 0, Rt = R%_i. 

Proof: With exactly active relays, and with k < n < N, the resulting rate can be written explicitly 
as: 

Rn = max min {(toLsr^), (toLsr2 + tlLrj^r2), ■ ■ ■ , (toLsr^ + ■ ■ ■ + tk-lLr^^irk): 
to,---,tN 

(toLsr„ + . . . + tk-lLr^_^rn + tkLr^r„ + tk+lLr^+irn ■ ■ ■ tn-lLr„_-^r„) ■, 
{toLsd + . . . + tk-lLr^_j^ra + tkLr^r^ + ^k+lLr^+ird + • • • ^Af-^rjvr^) } • 

Setting tk = gives 

Rn = max min {(to-^sn), (^0-^sr2 + ^l-^rir2)) • • • > (^O-^sr^ + • • • + i^fc-l-^rfe.irfc); 

tOv:tfe-l:tfe + l v-tjV 

(toLsrn + . . . + tk-lLr,._j^r„ + ^fc+l-^rfe+ir„ ■ • • ^n-l-^r„_ir„ ) , • • ■ , 
(toLsd + . . . + tk-lLr,,_^ri + tk+lLrf.^j^ra + • • • ^NLrj^rc)} (8) 
< max min {(to-^^^rj, (^O-^^sra + ^l-^^nra), • • • , 

to,...,tfc_i,tfc+l,...t]V 

itoLsr„ + . . . + tk^iLri^_j^r„ + ^fc+l-Z^rfc+ir„ • • • in-l-^r„_ir„ ) , • • • , 
(toLsd + . . . + tk-lLr^_^ri + tk+lLrf.^ir^ + ■ ■ ■ t^Lr^ra)} (9) 

= Rn-1, (10) 

since ^ has one fewer term in the minimization than ([8]). ■ 

To solve the optimization problem of (01) we thus require only the critical points for which tl ^ 0, V/c. 
In the following proposition, we show that for each Rj^j, i.e., given a set of potential relays, only one 
solution satisfies tl ^ 0,Wk. 

Proposition 2: The unique solution to the optimization problem in the minimization in (H]) for which 
tl ^ 0,\/k is given by Ji(ti) = hitiM) = ■■■ = /tvI^i, • • • ,t7v) = Initu.. . ,t7v)- 

Proof: We consider all possible critical points obtained from the optimization in (Hj). The points are 



. . . , 

(7) 
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obtained either by maximizing each individual term in (HI) or by intersecting all possible combinations of 
the terms in (Hj). We show that the only solution leading to non-zero solutions results from intersecting 
every term in (HI). 

The critical points for the optimization problem can be obtained by solving the following: 

1) Maximize the individual terms in (H)) except Id{to, . . . ,1^): 

\/k<N, max 4(to,...4-i) s.t. to + • • • + 4-i < 1- (H) 

Because the optimization is not over tm,Wk < m < N, the solution to this problem clearly has all 
tm = < m < N, and thus cannot be a solution to the overall optimization problem. 

2) Maximize /^(to, • • • , t^)'- 

max Id{to, . . . tiv) = max {to^sd + • • • + t^Lrj^ra), s.t. to + • • • + ^jv < 1- (12) 

<0,...,tjv to,...,tjv 

In this case, all variables are included in the optimization. It is easy to show, however, that this func- 
tion is maximized by selecting the largest L value, i.e., evaluating the Kuhn-Tucker conditions leads 
to a solution of the form tm = l,tk = 0, Wk ^ m, where m = arg maxk{Lsd, L^^d, • • • , -^r^d, • • • , Lrj^d} 
Therefore, this solution is also not a solution to the overall optimization problem. 

3) Maximize the function that results from the intersection of all possible combinations of the functions 
Ik- Let denote all possible subsets of {1 . . . A^}. ^A then contains 2^ such subsets, i.e., \ Ai \ = 2^. 
Consider one such subset Sk = (mi, m2, . . . , rrik), with mi < m2 < rrik. One critical point then is 

max /^^(to, • • -imfc.!) (13) 

such that 

^rrii (^0) • • • trrii-l) = Im2 {^0, ■ ■ ■ im2-l) = ■ ■ ■ = /m^ (^0; • • • ^mfc-l)- (14) 

This optimization then gets repeated for all sets 5k E J^.ln all but one combination, this optimization 
is not over all the variables {to, . . .In}- As in point (1), this maximization also leads to = for 
some value of k. 

4) Maximize the intersection of all terms in dH): 

h(to) = -^2(^0) tl) = ■ ■ ■ = ^Nito, ■ ■ ■ , tN-l) = Id(to, ■ ■ ■ , tN)- (15) 



This is the only case that leads to 7^ 0, Vfc = . . . A^. ■ 

Essentially, this proposition shows that if all relays are to contribute, all terms in the minimization 
in ^ must be equal. This proposition applies to any value of N. Therefore, if the optimal solution has 
k < N relays, an expression like (HJ can be written for those k relays. 



B. Optimal solution 

The linear system of equations in (fT5l) has a simple solution. Setting each equation to a constant, solving 
for the vector of unknowns t = (to . . .t^) and normalizing, we obtain 



LAf+1, 



|Ljv^^l7V+l||l 



1^ T -1 1 



(16) 



where ||v||i denotes the sum of the elements of v, i.e., the 1-norm. Iat+i is the length-(A^+ 1) vector of 
ones and Lat+i is the (A^ + 1) x (A^ + 1) rate matrix 
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Lr2r'i 
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Lsd 


Lrid 


Lr2d 


• Lr^d 



(17) 



The solution in (|T6l) does not guarantee that the constraint tk > Wk = . . . N is satisfied. To ensure 
that only solutions for which this constraint is satisfied are considered, we again consider the set Ai. 
Each entry in the set corresponds to a rate matrix, L^, similar to that in (fTTI) . formed using the relays 
in that entry of the set. Furthermore, let |m| denote the size of the rate matrix L^- A relay set and its 
corresponding solution, denoted as t^, is included as a potential solution if satisfies the constraint, 
i.e.. 



(18) 



where 0|m| is the all-zero vector of size |m|, 0|m| = [0,0,0, ...0]-'^ and the inequality operates on an 
element-by-element basis. Let the set JC form the subset of ^A that comprises all potential solutions. Let 
Lfe, tfc and |A;| denote the rate matrix, its corresponding solution and size, respectively, for each entry of 
the set /C. Note that the number of active relays being considered in each entry is \k\ - 1. Finally, the 



optimum solution can be obtained by solving (fT6l) for all possible combinations of active relays in the 
set /C i.e., 



t* = max ^'^.i'^' ,Vfc = 1,...,|/C|. 
Given that entry k* corresponds to t*, the maximum achievable rate vector can thus be written as 



(19) 



Li.«t* — Ilk* 



and the maximum achievable rate, R*, is 



l|fc*| 



(20) 



R* 



(21) 



Note that the solution described above is equivalent to the iterative maximization in and that removing 
a relay translates to removing the k^^ row and (k + iy^ column from the rate matrix in (fTTI) . Removing 
the first relay, for example, reduces the rate matrix in (fTTI) to 



■IN 







■'sd ^r2d 







Lr^d 



(22) 



Since \M \ = 2^, 2^ possible solutions must be tested to find the global optimum. 



C. Numbering 

In Section IIII-A[ we gave the solution to the optimization problem for a network with nodes numbered 
as in Figured! The numbering of the relay nodes impacts performance through causality: relay decodes 
information from relay rk-i, but not vice-versa. A complete solution to the optimization problem must 
therefore take into account an optimal numbering scheme. In the worst case (in terms of computational 
power), an optimal solution can be obtained for a specific numbering scheme, and the truly optimal 
solution can be maximized over all possible numbering schemes. 

Clearly, such an approach is impractical. Although a search for an optimal or effective sub-optimal 
solution is beyond the scope of this paper, we study the effects of numbering on the solution and resulting 
rate by considering some numbering schemes based on heuristics. We consider two approaches: numbering 
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based on average channel conditions, and numbering based on instantaneous channel conditions. 

1 ) Numbering based on average channel conditions: In the case of the linear network in Figure [B 
the numbering is trivial: node numbers increase away from the source and towards the destination. In 
the case of square network with nodes arranged in a grid, we consider two numberings which we refer 
to as Average Descending Numbering and Average Linear Numbering, shown in Figures [3] and [2l 
respectively, for a 4 x 4 network. 

• Average Descending numbering: node numbers increase towards the destination and downwards, 

• Average Linear numbering: node numbers increase towards the destination but vertical numbering 
ensures that nodes closest to each other retain close numbering. 

2) Numbering based on instantaneous channel conditions: 

• Instantaneous S — Rk numbering : node numbers increase with increasing source-relay channels. 
The first node has the best source-relay channel, the second node has the second-best source-relay 
channel, etc. 

• Instantaneous Rk — Rm numbering : nodes are numbered to maximize the channel between adjacent 
nodes. The first relay has the best source-relay channel. The second relay has the strongest ri -relay 
channel. Numbers are assigned in this process to unoccupied relays. This heuristic is based on the 
notion that we should maximize the capacity of each (R^, Rk+i) hop. 

• Random numbering : nodes are numbered randomly. This case evaluates the worst-case scenario 
and tests the robustness of the optimization to numbering. 

These schemes are evaluated via simulations in Section |Vl As we will see, the achievable rate is remarkably 
robust to the chosen numbering scheme. 

D. Partially Connected Network 

In this section we briefly discuss the more practical case of a partially connected network in which some 
links between the nodes in the network are unavailable. This is a generalization of the fully-connected 
network discussed in Section ITlI-AI above. Such a network is more likely to represent a large scale network 
where, in any case, the solution in (fT9l) would be computationally infeasible. 

As an example, consider the two-relay network with the link between ri and r2 is removed. The rate 
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matrix thus becomes 



-'sri 



L 








sr2 

Lsd Lr^ d L., 



(23) 



-'rid J-'r2d 

Removing the link thus reduces the rank of this matrix by one, and the rate matrix is now non-invertible, 
eliminating the solution defined by Ii = I2 = I3, where both relays are active. The optimal solution in this 
case is thus to select ri, r2, or not to relay. Note, however, that removing a link does not automatically 
lead to a non-invertible rate matrix. Consider, for example, the three-relay network with the link between 
ri and removed. The corresponding rate matrix 



Lsri 











Lsr2 


Lrir2 













Lr2rg 





Lsd 


Lrid 


Lr2d 


Lr^d 



(24) 



is full-rank and invertible. 

The approach to the optimization problem for the case of the arbitrary connected network is that the same 
as for the fully-connected network, with the exception that the rate matrix L^+i may not be invertible, 
in which case the corresponding solution is inadmissable. The remaining steps remain unchanged. 

IV. Implementation with Reduced Complexity 

The solution to the optimization problem in (HI), ([5]) and Q involves checking 2^ potential solutions. 
Although the process is conceptually simple, each solution involves the inverse of a rate matrix. In this 
section, we show how the optimization problem in the previous section can be significantly simplified 
using a recursive solution. This solution, which exploits the special structure of the rate matrix, greatly 
simplifies the matrix inversion, as well as reduces the number of possible solutions to check. Essentially, 
while the solution in Section ITlI-AI was a top-down approach, the approach we suggest here is bottom-up. 

Consider a set of p relays, V = {ri, r2, . . . , rp},p > 0, and its corresponding rate matrix lip^i, solution 
vector tp^^ and maximum rate (if available) . We note that if p = and the set is empty, the rate 
matrix and solution vector are constants, Lsd and 1, respectively. Denote as V the set V appended with 
another relay, i.e., V = {ri, r2, . . . , r^, rp+i}. Denote as Lp_j^2' ^^'^ LC^' the matrix, solution vector 
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and rate corresponding to set V'. 

Proposition 3: Given (Lp+i) ^ , (Lp+2) ^ can be obtained with computational complexity order of 

0{f) 



Proof: For p > 0, the rate matrix Lj^2 can be written 



as 



S+2 



Lj;i(l:p,l:p) 


Opx2 


F2XP 


'^2 



(25) 



where Lp+i(l : P, 1 : p) denotes the first p rows and columns of the rate matrix Lj+i, 0^x2 is a (p x 2) 
matrix of zeros, T2 is a (2 x 2) lower- triangular matrix, and F2xp is a (2 x p) fully-loaded matrix. Note 
that Lp+i(l : p, 1 : p) is triangular. Using the inverse of a partitioned matrix [39], (Lp+2) ^ can be written 
as 



(Lj;i(l:p,l:p)) 



-1 



-T2-'F2xp (Lj+,(l:p,l:p))~' 







px2 



rp-1 
^2 



(26) 



Note that (Lp^_]^(l : p, 1 : p)) ^ is the inverse of a partition of the triangular matrix Lp^_j^. Using the inverse 
of a partitioned matrix one more time, however, it is easy to see that 



and thus 



and hence obtaining (Lp+2) ^ OijP') operation. ■ 

Using this above proposition, the solution vector t^^2 ^^'+2 can be obtained from the solution vector 



. :p, 1 


:p))" = (L?-m) 


"(1: 


p, 1 : 






(L?+i)-^(l:P,l 


: p) 




0px2 


-T2- 


^F2xp(Lj+J-Hl 


:p, 1 


:p) 


-•-2 



(27) 



(28) 



p+2 



-^p+2 



(Lp|2) lp+2 



t?+i(l:p) 



t?;2(p+i) 

Ltj;2(pH-2) 



(29) 



where t?^i(l : p) represent the first p entries of the already-calculated solution vector t^^^, and t?|.2(p+ 1) 
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and tp_j_2(p + 2) are the last two entries of the solution vector tj^g that remain to be calculated. R^' — 
is the maximum achievable rate obtained using the set V' of relays. The last two entries 



of the solution vector tp^2{P + 1) ^p+2{P + 2) can be written as 



R 



t?;2(P+2) 

With a corresponding achievable rate R^' 

1 



T-^Fsxp (Lj+J ' (1 : p, 1 : p) | T^-^ ] l(p+2)xi, 



(30) 



R^' = 



(31) 



where we use j) to denote the summation over all the elements of matrix A. 

Using the above, the optimization problem for a network of N potential relays can be solved recursively 
as follows: 

1) Determine the set of all potential relay combinations. Sequence the set as: 

M = {(ri), (ri, ra), (n, ra, rg), . . . (n, rs, . . . , rjv), 
(ri, ra), (n, rg, r4), . . . , (n, rg, . . . , tat), 



(n,rjv), 



(r2), (r2, ra), (ra, rg, r4), . . . (rs, rg, . . . , tat), 
(r2, r4), (r2, r4, rg), . . . , (ra, r4, . . . , tat), 



('^2,rAr), 



(rAr-i,riv)}. 
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Note that each "row" of is a subset of relay combinations in which each element is formed from 
the previous element by adding a relay. 

2) In each "row", obtain the rate matrix, its respective optimized time allocation vector and achievable 
rate for each element (i.e., relay combination) recursively using (|28] ). (|29l ). (|30l ) and (|3T1) . 

3) Check that for each particular set V of p relays, the solution tp and achievable rate Rp satisfies the 
constraints: 

R'P > 0, (32) 
tJVi > Op+1. (33) 

• If both constraints are satisfied, place the solution in the potential set of valid solutions /C, 
advance elements and return to step (1). 

• If (|33l) is not satisfied, check which element of the the allocation vector tp does not satisfy the 
constraint. 

- If any of the first (p — 1) entries of tp are less than zero, i.e., tp(l : p — 1 < Op_i), this 
constraint will not be satisfied for any other relay combinations in this "row". Advance rows 
and return to item (1). 

- If the constraint is not satisfied by either of the last two items in the solution vector, discard 
the solution but check the other elements in the "row". 

4) From the set /C, pick the highest achievable rate and its corresponding time allocation. 
The recursive algorithm given above simplifies the optimization problem in two ways: 

1) It reduces the computation load of determining successive matrix inverses by writing each matrix 
inverse as a function of another, already known, matrix inverse, and two other matrices obtained 
through simple matrix multiplication. 

2) It may eliminate infeasible solutions by discarding relay combinations which do not satisfy con- 
straints. For example, if the relay combination (ri,r2,r3) does not satisfy the constraints, the 
combination (ri, r2, rs, r4) may be automatically discarded. 

A. Complexity and Number of Operations 

In the next paper we will compare relay selection schemes partly on computational complexity. In this 
section we calculate this complexity, which also quantifies the computational savings of the recursive 
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scheme presented above in Section |Wl 

The complexity of the recursive scheme is bounded by complexity of matrix multiplication. The number 
of operations (multiplications and additions) required in the product of two matrices of size (m, n) and 
(n, p) is 2mpn — mp [40], and the number of operations required for the product of a matrix of size 
(m, n) with a square, size-n diagonal matrix is 

^n— 1 n 



m 



^fc + ^fc = mnl (34) 

,fc=0 k=l / 

We now calculate the number of operations required for each rate matrix of size (g + 1), corresponding 
to the set Q' of q relays. The calculation of the matrix fundamental to the recursive algorithm, 

(lSi)"' = -T^ iF2x(,-i) {Lfy' (1 : g - 1, 1 : g - 1). (35) 

requires a total of + 2q + 1 operations, broken down as: 

1) — T2"^ 5 operations, 

2) — T2""'^F2x(q-i) = A2x(g-i) ^ 6(g — 1) operations using 2mpn — mp, 

3) A2x(g~i) (L^)"^ (1 : g - 1, 1 : g - 1) ^ 2(g - 1)2 operations, using dMl). 

From (|3T1) . the number of operations required to calculate R^' is g^ + 2g + 4. Using (l30l) . the number of 
operations required to update the solution vector is 1 + 2(g + 1) = 2g + 3. Summing the above, we obtain 
the total number of operations required in one iteration of the resource allocation algorithm: 

Op(g) = (2g2 + 2g + 1) + (g^ + 2g + 4) + (2g + 3) = 3q^ + 6g + 8. (36) 

Note that the complexity order of calculating each rate and solution vector is 0{q^). Without the 
recursion, this complexity is of order 0{q^), resulting from the inverse of the rate matrix. The recursion 
thus introduces significant savings in terms of complexity. 

We now calculate the worst-case total number of operations required by the resource allocation algo- 
rithm. In the worst case, the algorithm cycles through 2^ operations consisting of (^) sets of g relays 
which require 3q'^ + 6g + 8 operations. The total worst case number of operations is therefore 

J(3g2 + 6g + 8). (37) 

This calculation could be rendered more precise if it were possible to account for the savings obtained in 
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Section |IV] which eliminates some infeasible solutions a priori by discarding relay combinations known 
to not satisfy the constraints. The probability of this occurring for particular channel realizations is 
unfortunately very difficult to compute, and we thus show only the worst-case result. 

V. Simulations 

In this section, we present results of the resource allocation scheme discussed in Section [ni] for networks 
with 1 to 6 relays arranged linearly, and 4 and 9 nodes arranged in a grid. The figure of merit is the 
achievable rate Ra with an outage probability of 10^^, i.e. Fi[R* < Ra] = 10^^. A closed form expression 
for the outage probability of optimized cooperation is very complicated and beyond the scope of the paper. 
The outage probability and rate are thus obtained numerically. 

The relays are equispaced on a line between the source and destination, as in Figure [H and we use an 
attenuation exponent of = 2.5. This choice is motivated by the application of static mesh-nodes installed 
on posts; transmissions between such nodes should undergo little shadowing and a lower attenuation 
exponent. From 60000 fading realizations we obtain the cumulative density function of the instantaneous 
rate Fji(r). The outage rate is the rate for which the probability of outage is 10^^, i.e., F^^(10^^). 

Figure |4] and |5] plot the outage rate as a function of the average end-to-end SNR, for optimized 

and non-optimized cooperation, respectively. The rate for the optimized cooperation is obtained from (fT9l) . 
Non-optimized cooperation uses equal time allocation, i.e., the rate for a particular relay set is simply the 
minimum of the mutual information at each node. Non-optimized cooperation, however, does optimally 
select relays by choosing the best, in terms of rate, of the 2^ relay combinations. Comparing Figure |4] 
and Figure [5] shows that optimizing resources increases rates significantly, as expected. The outage rate 
increases as a function of nodes available to relay. We also note the typical phenomenon of decreasing 
marginal returns: the gains of adding each additional relay decreases with increasing number of relays. 

Figures [6] and |7] show the average number of relays that are active from the set of potential relays for 
optimized and non-optimized cooperation. For each network size, this number is a decreasing function. In- 
terestingly, the number of active relays decreases much faster for non-optimized as compared to optimized 
cooperation, suggesting that optimizing resources distributes the relaying burden more effectively. 

To test the effect of geometry on the outage rate, we compare the rates obtained by optimizing resources 
and the placing relays on a line, as in Figure |4] to those obtained by placing the relays on a regular square 
grid. We number the relays in the grid in ascending order downwards and towards the source; a derivation 
of the optimal numbering is beyond the scope of this paper. The results are demonstrated in Figure [H 
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where we place 4 and 9 relays on a 2 x 2 and 3x3 square grid. As shown in the figure, the rate for the 
linear constellation is significantly higher than that obtained by the grid constellation, suggesting that the 
path-loss incurred by traversing all the nodes laterally results in non-negligible performance loss. 

We evaluate the performance of the numbering schemes discussed in Section IIII-CI in Figure [9l The 
four schemes, including two based on average channel conditions and two based on instantaneous channel 
conditions, exhibit indistinguishable performance in terms of rate. There is an expected drop in rate with 
random numbering, though we note that this drop is no more than approximately 0.25 bits/channel use. 
The algorithm is thus quite robust to numbering schemes. 

Figure |9] also shows the outage rate for a network with randomly placed nodes. Here the node locations 
are chosen from a uniform distribution over an area equivalent to that of the square grid . The internode 
channels are obtained as before. This example eliminates possible dependencies of the results obtained 
earlier on the chosen array geometry. The numbering here is based on the instantaneous S — Rk channels. 
In such a random network, as expected, the available outage rate is lower than in a square grid network; 
however, at higher SNR levels this difference disappears. Again, the significant gains due to resource 
allocation are clear. 

In Figure [9] we also compare the effect of numbering when used without resource allocation, and 
show only the case of instantaneous S — Rk numbering and random numbering. The improvement from 
instantaneous over random numbering in this case is less than 0.1 bits/channel use. The robustness of the 
numbering scheme thus increases by eliminating time optimization. To gain insight into this phenomenon, 
in Figure \T0\ we plot the average number of active users for the instantaneous and random numbering 
schemes with and without resource allocation. We first observe that the instantaneous numbering scheme 
uses more relays than the random numbering scheme when resource allocation is used, and that this 
difference is constant over the SNR region of interest. Without resource allocation, on the other hand, the 
number of relays used when using instantaneous and random numbering decreases quickly and is constant 
for SNR values higher than 10 dB. It is clear from this figure that the difference in rate performance 
between instantaneous and random numbering is an increasing function of the number of selected relays. 
Because so few relays are selected without resource allocation, the effect of the numbering scheme is 
negligible. The influence of the numbering scheme increases when time allocation is introduced, increasing 
the number of relays used for both numbering schemes and increasing the sensitivity to the numbering 
scheme. This sensitivity increases slowly, however, and is negligible for the various numbering schemes 
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based on heuristics. 

VI. Conclusions 

In this paper, we determined the optimal channel resource allocation, in terms of time allocation, for the 
N-node cooperative diversity, multihop network using DF and independent codebooks. For a particular 
network, i.e., set of potential relays, the unique solution for a particular relay numbering scheme is obtained 
by taking the inverse of the triangular rate matrix, and the optimal solution is found by maximizing 
over the rate for each possible network, given its maximum size. Through simulations, however, the 
optimization is shown to be robust to the numbering scheme. We show that by exploiting the special 
structure of the rate matrix, the optimization can be performed in a recursive fashion which decreases 
the computation load of the rate matrix inverse and the number of required iterations. Node selection is 
inherent to the optimization strategy. Simulation results show significant gains in achievable rate due to 
resource allocation, but diminishing marginal returns as a function of network size. Furthermore, we show 
a significant benefit to arranging the nodes in a linear, as opposed to a grid, constellation. 
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Fig. 1. Location of the relays with respect to the source and destination. 
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Fig. 2. Numbering in a square 4x4 network in linear order 
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Fig. 3. Numbering in a square 4x4 network in descending order 
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Fig. 4. Outage rate vs. SNR using 1, ... 6 potential relays and with resource allocation. 
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Outage rate with no resource allocation for 1 to 6 relays arranged in a line 
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Fig. 5. Outage rate vs. SNR using 1, ... 6 potential relays and without resource allocation. 
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Average number of active relays with resource allocation for 1 to 6 Relays arranged in a line 
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Fig. 6. Average number of active relays with 1, ... 6 potential relays and with resource allocation. 
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Fig. 7. Average number of active relays with 1, ... 6 potential relays and without resource allocation. 



Outage rate with resource allocation for 4 and 9 relays arranged in a grid and in a line 
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Fig. 8. Outage rate vs. SNR using resource allocation and for 4 and 9 relays arranged in a grid and in a line. 



Outage rate w ith 9 nodes 

6 Square Network. Instantaneous S-R^ 

Square Network. Instantaneous Rk-Rm 
Square Network. Average Linear 




-10 -5 5 10 15 20 

SNR (dB) 



Fig. 9. Outage rate vs. SNR using resource allocation and for various numbering schemes for 9 potential nodes arranged in a grid 
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Fig. 10. Average number of active relays using resource allocation and with various numbering schemes for 9 potential nodes arranged in 
a grid 



